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ABSTRACT 

We report precise Doppler measurements of the nearby (d — 10.34 pc) M dwarf G1649 that reveal 
the presence of a planet with a minimum mass Mp sin i = 0.328 Mj up in an eccentric (e = 0.30), 
598.3 day orbit. Our photometric monitoring reveals G1649 to be a new variable star with brightness 
changes on both rotational and decadal timescales. However, neither of these timescales are consistent 
with the 600-day Doppler signal and so provide strong support for planetary reflex motion as the best 
interpretation of the observed radial velocity variations. Gl 649 b is only the seventh Doppler-detected 
giant planet around an M dwarf. The properties of the planet and host-star therefore contribute 
significant information to our knowledge of planet formation around low-mass stars. We revise and 
refine the occurrence rate of giant planets around M dwarfs based on the California Planet Survey 
sample of low-mass stars (M* < 0.6 M Q ). We find that / = 3A±l^% of stars with M* < 0.6 M© 
harbor planets with Mp s'mi > 0.3 Mj up and a < 2.5 AU. When we restrict our analysis to metal-rich 

stars with [Fe/H] > +0.2 we find the occurrence rate is 10.7^42%- 

Subject headings: techniques: radial velocities — planetary systems: formation — stars: individual 
(G1649) 



1. INTRODUCTION 



Compared to the knowledge gleaned from the large 
sample of giant planets around Sun-like stars, little is 
known about the characteristics of Jovian planets around 
M dwarfs. This is due primarily to the empirical find- 
ing that the occur rence rate of detectab le planets scales 
with stellar mass ([Johnson et al.ll2007t ): low-mass stars 
(M* < 0.6 Mfn) simply do n ot harbor giant pl anets very 
frequently (jEndl et all2003l:lButler et alJ2006f >. The fre- 
quency of giant planets with Mp sin i > 0.3 Mj UD around 
Sun-l ike stars is 8% within 2.5 AU (|Cumming et al.1 
2008), and the occurrence of giant planets around M 
dwarfs is roughly a factor of 4 lower. 

While the lower masses of M dwarfs decreases the 
likelihood of giant planet occurrence, a handful of Jo- 
vian planets have been discovered around low-mass stars. 
The sample of M dwarfs known to harbor at least one 
Doppler-detected giant planet (Mp sin i > 0.2 Mj up ) 
is listed in Table [T] and shown in the H-R diagram in 
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Fig. 1. — Our Keck sample of low-mass stars plotted in the {V — 
K, Mk} plane. The solid line is a fifth-order polynomial fit to 
the mean main sequence for stars within 10 pc, which Johnson 
& Apps (2009) identify as an isometallicity contoure with [Fe/H] 
equal to the mean value of the Solar neighborhood. The dashed line 
corresponds to [Fe/H] = +0.2 based on the calibration of Johnson 
& Apps. The five-point stars show the positions of all of the M 
dwarfs known to harbor at least one giant planet. The solid circle 
denotes the position of Gl 649. 



Figure [TJ Also given in that table are the stellar and 
planetary masses from the literature, and stellar metal- 
licities from the broad- band photometric calibration of 
Uohnson k. Appd (|2009l ). 

These planets and their host stars demonstrate that 
stellar mass is not the only characteristic that corre- 
lates with the probability of a star harboring a planet. 
Stellar metallicity has been shown to be a strong 
predictor of planet occu rrence around Sun-like stars 
(|Fischer &: Valentlll2005af ). and the correlation between 
planet frequency and stellar m etal content app e ars to 
hold for the M dwarfs, as well. Uohnson fc Appsl (|2009f ) 
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Fig. 2. — Our Keck/HIRES radial velocity time series for G1649. 
scatter of the residuals (bottom panel) is 4.2 m s — 1 , and yxF = 1.17. 

found that M dwarfs with Jovian planets tend to be sig- 
nificantly metal-rich compared to a 10 pc, volume-limited 
sample of stars on the lower main sequence. For example, 
Gl 849 harbors a Jovian planet in a long-period orbit and 
is among the most metal-rich stars in the Solar neighbor- 
hood with [Fe/H] > +0.45. 

If this preliminary trend proves to be real then it 
will provide valuable constraints for theoretical mod- 
els of planet formation around a broad range of stel- 
lar characteristics. The effect of metallicity on planet 
occurrence will also inform the target selection of fu- 
ture Doppler and tran sit surveys targeting low-mass stars 
(e.g. llrwin et al.f2008l ). as well as the interpretation of re- 
sults from dir ect-imaging, astrometr i c and microlensing 
surveys (e.g. INielsen fc Close! l20"09t IPravdo fc Shaklanl 
120091: bong et al.ll2009D . ~ 

As the time baselines, sample sizes and Doppler preci- 
sion increase for the various Doppler surveys of low-mass 
stars, the relationships between the physical character- 
istics of stars and the properties of their planets will 
come into sharper focus. We are monitoring a sample 
of 147 late K and early M stars as part of the Cali- 
fornia Planet Survey at Keck Observatory with a cur- 
rent temporal baseline of « 12 years and Doppler pre- 
cision of 2-3 m s _1 (Johnson et al. 2007; Howard et 
al. 2009b). In this contribution we announce the detec- 
tion of a new Saturn-mass planet orbiting a nearby M 
dwarf. Gl 649 is only the eleventh M-type star known to 
harbor at least one Doppler-detected planet 10 , and it is 

10 Several additional low— mass host stars have been discov- 
ered by gravitational microlensing surveys (Bond ct al] 120041 ; 



rhe dashed line shows the best-fitting Keplerian model. The rms 

only the seventh low- mass star with a Doppler-detected 
giant planet (see also IBonfils et al.l 120051: IManess et al.l 
I2007t iForveille et al.|[2009t iMavor et alJl2009L for exam- 
ples of low-mass planet detections) . In the following sec- 
tion we describe the stellar properties of G1649, and our 
spectroscopic observations and Doppler-shift measure- 
ments. In §[3] we test the validity of our interpretation of 
the observed radial velocity variations by measuring the 
false-alarm probability and by examining our photomet- 
ric measurements. We conclude in § f4] with a summary 
and discussion of G1649b, and we place this latest ex- 
oplanet in context with other giant planets discovered 
around M dwarfs. 

2. OBSERVATIONS AND ANALYSIS 

2.1. Stellar Properties 

G1649 (= HIP 83043) is an Ml. 5 dwarf with a Hip- 
parco s parallax-ba s ed di stance of 10.34 ± 0.15 par- 
sees (jvan Leeuwenl 120071) . apparent magnitude V — 
9.7, a nd absolute magnitude My = 9.627 ± 0.053 
(|ESAI Il997t) . We use the b r oadba nd metallicity 
calibration of Uohnson fc Appsl ([2009) to estimate 
[Fe/H] =+0.08±0.06, an d we adopt the stella r mass es- 
timate provided by the iDelfosse et al.l ((2000) K s -band 
mass-luminosity relationship, which gives M+ = 0.54 ± 
0.05 M w . Using the infrared flux method. lAlonso et al.l 
( 1996 ) give an effec tive temperature T g = 3700 ± 60 K. 
Wrig htet al.l (|2004f ) measured the emission in the Call H 
emission line rela tive to the stellar p hotosphere on the 
Mt. Wilson scale {D uncan et al.lll991J ) and give a median 

IGould et ani200"l [Beaulieu et al.ll200"l IDone et al.ll2009l. e.g.). 
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Fig. 3. — Periodogram analysis of our RV time series, which 
reveals a strong peak at P = 592 days. The dashed lines show 
the analytic false-alarm probability of a peak arrising from noise 
sampled at our times of observation. 

"grand S" value of 1.55. This S value places the chromo- 
sphcric activity of G1649 among the top 20% of nearby 
early M-type stars, as shown by iRauscher fc Marcvl 
(2006) and by lGizis et all (|2002l ). Our spectra show H- 
alpha to be in absorption, as was found for all B aimer 
lines observed in the spectrum of G1649 jG izis et all 
2002) . The stellar properties of Gl 649 are summarized 
in Tabled 

2.2. Radial Velocities and Keplerian Fit 

We began monitoring G1649 at Keck Observatory in 
1999 Oct ober using the H igh-Resolution Echelle spec- 
trometer (|Vogt et al.l [19941 in our standard Iodine cell 
setup with the B5 decker, giving a reciprocal reso- 
lution A/AA = 55, 000 p er ~ 4-pixel resolution ele- 
ment ([Howard et al.ll2009f ). We measured the Doppler 
shifts of the star from each star-times-iodine observa- 
tion using the stan dard analysis procedure presented by 
iButler et al.l (|T996), with subsequent improvements over 
the years. For HIRES observations made prior to the 
2004 CCD upgrade the measurement uncertainties range 
from 3.3-4.4 m s , and improve to 0.9-1.5 m s _1 there- 
after. 

Our 44 radial velocities are presented in Table [3] (with- 
out jitter) and the time series is shown in Figure [5] 
(with jitter). The scatter in the measurements is larger 
than expected from the measurement errors, and a peri- 
odogram analysis of the data reveals strong power at pe- 
riods near 592 days, with a corresponding analytic false- 
alarm probability < 0.0001 (Figure 03). 

We used the partially -linearized Keplerian fi tting code 
RVLIN 11 described by I Wright fc Howard! (|2009[ ) to search 
for a best-fitting orbital solution to the data. To en- 
sure proper weighting of our measurements in the fitting 
procedure, we inflated the error bars to account for RV 
noise from astrophysical sources. This stellar "jitter" 
term is calculated based on the star's chromospheric ac- 
tivity, B — V c olor and absol ute V-band magnitude using 
the formula of lWrightj |2005). We adopt a jitter estimate 

11 http://cxoplancts.org/codc/ 



of 3 m s _1 for G1649, which we add in quadrature to the 
measurement errors. 

We find that a single-planet Keplerian model with a 
period P = 598.3 ± 4.2 days, eccentricity e = 0.30 ± 0.08 
and velocity semiamplitude K = 12.4il.lms" 1 results 
in a root-mean-squared (rms) scatter of 4.2 m s _1 in 
the residuals and \/aJ = 1.17, indicating an acceptable 
fit 12 . The resulting minimum planet mass is Mpsini = 
0.328 Mj up , and the semimajor axis is a — 1.135 AU. 
The best-fitting solution is shown in Figure [2l with the 
residuals to the fit shown in the lower panel. The orbital 
parameters are listed in Table [D 

The parameter uncertainties given above were esti- 
mated using a Markov Chain Monte Carlo (MCMC) 
algorithm with 10 7 links, in which a single randomly- 
chosen parameter was perturbed at each link, with a per- 
turbation size tun ed such that 20-40% of the ju mps were 
executed (see e.g. lR^l27MrWinn et al.H20"08l and ref- 
erences therein). The resulting "chains" of parameters 
form the posterior probability distribution, from which 
we select the 15.9 and 84.2 percentile levels in the cumu- 
lative distributions (CDF) as the "one-sigma" confidence 
limits. In most cases the posterior probability distribu- 
tions were approximately Gaussian. 

3. TESTING THE NULL-HYPOTHESIS 

3.1. False-Alarm Probability 

The Doppler semi-amplitude of our best-fitting model, 
K = 12 A m s j is comparable to the measurement un- 
certainties and stellar jitter, which prompted us to test 
the null-hypothesis that the apparent periodicity arose 
by chance from larger-than-expected radial velocity fluc- 
tuations and sparse sampling. We tested this possibility 
calculating the false-alarm probability (FAP) based on 
the goodness of fit statistic Ay 2 ([Howard et al.l 120091 : 
IMarcv et all 120051: [Cummind 120041) . which is the differ- 
ence between two values of xt'- one from the single-planet 
Keplerian fit and one from the fit of a linear trend to the 
data. Each trial is constructed by keeping the times of 
observation fixed and scrambling the measurements, with 
replacement. We record the Ay 2 value after each trial 
and repeat this process for 10,000 trial data sets. For 
the ensemble set we compare the resulting distribution 
of Ay 2 to the value from the fit to the original data. 

We found that none of the 10 4 trials resulted in a higher 
value of Ay 2 ,, which we interpret as a < 0.0001 probabil- 
ity that the 600-day periodicity is a spurious signal due 
to random fluctuations. 

3.2. Photometric Variability 

We note that our FAP value only addresses the exis- 
tence of a periodicity in the radial velocities, but does 
not test its cause. As an additional test of the null hy- 
pothesis we acquired brightness measurements of Gl 649 
in the Johnson V passband with the T3 0.4 m auto- 
matic photometric telescope (APT) at Fairborn Obser- 
vatory. The APT observations cover five observing sea- 
sons between 2001 February and 2009 June and reveal 

12 We use \/x5 to indicate the factor by which the observed 
scatter about the best-fitting model differs from our expectation 
based on the measurement errors. Thus, the scatter about our 
model is a factor of 1.17 larger than our average error bar. 
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Fig. 4. — Top Panel: Our 337 V — C photometric observations 
of Gl 649 in the Johnson V band acquired with the T3 0.4m APT 
at Fairborn Observatory. Second Panel: The K — C observa- 
tions plotted with an identical scale as the top panel. Comparison 
of these two data sets shows that Gl 649 varies in brightness on 
night-to-night and year-to-year timescales. Third Panel: Obser- 
vations from the second observing season plotted on an expanded 
x axis clearly show low-amplitude brightness variability in Gl 649. 
Fourth Panel: Frequency spectrum of the observations from sea- 
son two gives a best period of 24.55 days. Bottom Panel: Plot 
of the data from season two phased with the 24.55-day period re- 
veals coherent variability with a peak-to-peak brightness amplitude 
0.012 mag. 



photometric variability in Gl 649 on both rotational and 
decadal timescales. Details on APT operations, data 
acquisition and reduction proce dures, and p r ecision of 
the observa t ions c a n be found inlHenrv et all (|1995bl lal); 
iFekel et al.1 (|2005h : lEaton et alJ (|2003D . 

Our 337 Gl 649-minus-comparison (V — C) differential 
magnitudes are plotted against heliocentric Julian Date 
in the top panel of Figure 0J The comparison star is 
HD 152342 (V = 7.10, B-V = 0.35, F2V). Most obvi- 
ous in this plot are the year-to-year changes in the mean 
magnitude of the observations. The mean magnitudes 
have a range of 0.0126 mag and suggest the possible ex- 
istence of a spot (magnetic) cycle in Gl 649 with a length 
of at least several years (see, e.g.. lHenryl ll999: Hall et ahl 
2009). The top panel also shows that the range in the 
V — C observations is ~0.02 mag within all five observing 
seasons. The standard deviations of the individual five 
seasons are all between 0.0057 mag and 0.0063 mag. 

The check minus comparison (K — C) differential mag- 



nitudes are plotted in the second panel of Figure 2] at the 
same magnitude scale as the top panel. The check star is 
HD 153897 (V = 6.57, B-V = 0.43, F4V). The yearly 
means of the K — C measurements have a range of only 
0.0016 mag; their standard deviation is only 0.0007 mag. 
The standard deviations of the individual five seasons 
are all between 0.0034 mag and 0.0042 mag. Thus, the 
larger night-to-night scatter of the V — C measurements 
and the observed year-to-year change in the V — C mean 
magnitudes must both be intrinsic to G1649. 

The observations from season 2 are replotted in the 
third panel of Figure IH again on the same magnitude 
scale as in panels 1 and 2. Low-amplitude variability is 
clearly seen with a period of 20-30 days. Panel 4 shows 
our computed frequency spectrum for season 2, where 
the y axis plots the reduction factor in the variance o f 
the observations for each trial frequency (jVam cck 1971). 
We find a clear period of 24.55 days; the observations are 
plotted phased with this period in the bottom panel of 
Figure [4] A least-squares sine fit to the phase curve gives 
a peak-to-peak amplitude of 0.01 17±0. 0015 mag. Similar 
analyses yield periods of 23.72, 27.92, 25.80, and 21.86 
days for seasons one, three, four, and five, respectively. 
The mean of these five periods is 24.8 ± 1.0 days, which 
we take to be the rotation period of Gl 649 revealed by 
rotational modulation in the visibility of cool starspots 
on the photosphere of Gl 649. This photometric variabil- 
ity is consistent with the level of chromospheric activity 
i n the star, as mention ed above. 

iHenrv et al.l (|1995bf ) show many examples of ac- 
tive stars with low-a mplitude starspot variability. 
iQueloz et ail (|2001l) and lPaulson et alJ (2004(1 show sev- 
eral examples of stellar spots masquerading as planets. 
In the case of Gl 649 described here, the photometric ob- 
servations reveal variability timescales that are inconsis- 
tent with the 600-day radial velocity variations. Sinilarly, 
we examined the time- variability chromospheric emission 
from each of our spectroscopic observations. While the 
S- value, as measured from the Call H&K emission, has 
a variance of 0.14 dex, we observed no peridicities near 
the putative orbital period. The lack of photometric and 
chromospheric variability provide additional strong sup- 
port for the interpretation of planetary reflex motion as 
the cause of the observed radial velocity variability in 
G1649. 

4. DISCUSSION 

We have presented the discovery of a Saturn-mass 
planet (Mp sin i = 0.328 Mj up ) orbiting the nearby, 
low-mass star G1649 (d = 10.34 pc, 0.54 M ). G1649b 
resides in an eccentric (e = 0.30) orbit with a pe- 
riod of 598.3 days, corresponding to a semimajor axis 
a = 1.135AU. 

G1649 is only the seventh M dwarf known with a 
Doppler-detected giant planet, and the fifth detection 
from among the 147 low-mass stars we've monitored over 
the p ast decade at Keck Observatory (e.g. I Johnson et aD 
I2007f ). The low Doppler amplitude of the planet (K = 
12.4 m s _1 ) highlights our need to attain high measure- 
ment precision to find low-mass planets, and to maintain 
that precision over long time baselines to detect planets 
at larger semima j or ax es. 

John son et al.l (|2007f) recently analyzed the detection 
rate among our sample of low-mass stars and reported a 
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1.8% occurrence rate of planets with a < 2.5 AU. The 
increased time baseline and new detections of our sam- 
ple suggest that a reanalysis of the frequency of plan- 
ets around M dwarfs is warranted. Following John- 
son et al. we first note that the 10 year time 
baseline of our survey, together with our radial ve- 
locity precision, provides us with sensitivity to plan- 
ets with Mp sin i > 0.3 Mj up out to semimajor axes 
a w 2.5 AU, assuming and average stellar mass = 
0.5 M©. Note that in the analysis that follows, we ex- 
clude the recently detected planets G1832b, which was 
discove red by the Anglo-A ustralian Observatory planet 
search (|Bailev et al J 120081 ). and HIP 79431b, which was 
only recently added to the Keck survey as part of the 
metallicity-biased M-to-K program (Apps et al. 2010, 
submitted) . 

The probability density function (PDF) for the frac- 
tion of stars with planets given our number of detections 
k = 5 and total sample size N — 147 is given by the bino- 
mial distribution P(/|fc, TV) oc / 5 (1-/) 147 - 5 . Theoverall 
occurrence rate from our sample is given by the maxi- 
mum of the PDF, which we measure to be / — 3.4^'!%, 
where the upper and lower limits represent the 68.3% 
confidence interval measured from the cumulative distri- 
bution function. 

The corresponding giant planet fract ion around Sun- 
like s tars was recently measured by ICumming et al.l 
(2008, cf their Table 1), who report / = 7.6 ± 1.3%. In a 
similar study Bowler et al. (2009, submitted) measured 
the planet fraction around stars with M+ > 1.5 M© to be 
26^g%, albeit for minimum masses Mp sin i > 1 Mj up . 
Thus, the detection rate of giant planets around M 
dwarfs consistently lags behind that of higher mass stars, 
despite the enhanced detectability of planets around less 

massive stars since K cx M+ 2 ^ 3 for a fixed planet mass 
and period. The contrast between the measured planet 
fractions between M dwarfs and massive stars points 
to an even stronger correlation betwee n stellar mass 
and p lanet occurrence than measured by Uohnson et al.1 
(120071) . 

The correlation between stellar mass and planet for- 
mation is an important piece of observational evidence 
in support of the core accretion model of planet forma- 
tion. In this model, giant planets form in a bottom- 
up process, starting with the collisions of small dust 
grains and proceeding up through the formation of large 
proto planetary cores (see llda fc Linll200l lAlibert et all 
2005, for reviews). Once these cores attain a critical 
mass of ~ 10 M©, they can rapidly accrete gas from 
the surrounding disk. Given the limited lifetime of the 
gas disk, which dissipat es on timescales less than 5 Myr 
(jHernandez et al l 120081 : iCurrie et al1l2009D . the forma- 
tion of gas giant planets is a race against time that is 
rarely won in the protoplanetary disks of low-mass stars. 
The low density of raw materials, low orbital frequencies 

1/2 

(Q oc 1/P oc MJ at fixed a), and unfavorable temper- 
ature profiles in the disks around M-type stars greatly 
inhibit the core growth, which results i n a lower occur 
rence of giant planets ( Laughli n et al.l l2004t llda fc Lin 



2005 ; Kennedy fc Kenvonl 120081 : iDodson- Robinso n et al 



2009). 



Another importan t predictor of planet occu rrence is 
stellar metallicity. iFischer fc Valenti (|2005bf ) showed 



the fraction of Sun-like stars with planets correlates 
strongly with [Fe/H] , with an occurrence rate of ~ 3% for 
[Fe/H] < and a rise to » 25% for [Fe/H] > +0.3. Until 
recently, it was difficult to properly account for metallic- 
ity among the M dwarfs because the LTE spectral anal- 
ysis tools used for more massive stars are not amenable 
to th e complex spectra of low-mass stars ()Maness et al.l 
|2007|) . Because of the lack of knowledge about the metal- 
licity distribution of M dwarfs in general, and low-mass 
stars with planets in particular, it was difficult to deter- 
mine whether stellar mass or metallicitiy lay at the root 
cause of the puacity of planets around M dwarfs. 

The mass/met a llicity issue was recently addressed by 
Uohnson fc Appsl (|2009l ), who derived a revised broad- 
band photometric metallicity calibration for M dwarfs. 
They examined a sample of M dwarfs with F, G and K 
wide binary companions. By anchoring the metallicity 
of the M dwarf to its earlier-type companion, Johnson 
& Apps observed that metal-rich M stars reside "above" 
the mean main sequence of the solar neighborhood when 
viewed in the {V — K, Mk) plane. Further, they noticed 
the majority of the 7 planetary systems (containing plan- 
ets of all masses) that were known at the time contain 
metal-rich host stars. 

Figure [T] shows that with the addition of 3 new planet- 
host stars since the study of Johnson & Apps the planet- 
metallicity correlation among M dwarfs appears to per- 
sist. We can quantify this relationship by examining the 
fraction of stars in our Keck survey with [Fe/H] > 
that harbor giant planets. Including G1649, we find 
that all 4 of the stars that harbor at least one giant 
planet 13 fall within the subsamlple of 80 targets with 
[Fe/H] > 0. Based on this, we measure a planet fraction 
/ = 5.5^2 i% f° r [Fe/H] > 0. If we restrict our analysis 
to [Fe/H] > +0.2 (dashed line in Figure [J), 3 of these 33 
"super-metal-rich" stars harbor planets, corresponding 
to / = 10.7±||%. 

The uncertainties in our measured planet fractions are 
large due to the small sample sizes involved. This un- 
derscores the need for extending the time baseline of our 
current survey and expanding the target list to include 
additional low-mass stars. Future surveys of nearby, 
low-mass stars such as the M2K planet search (Apps 
et al. 2010, subm itted) and the MEarth transit survey 
(|Irwin et al.1 120091 ) will build upon the current sample 
and provide a clearer picture of the planet-metallicity 
relationship suggested from our analysis. A larger sam- 
ple of planets detected around M dwarfs will also pro- 
vide crucial leverage in understanding the relationship 
between stellar mass and planet properties, especially 
when compared to the growing sample of planets dis- 
covered around massive stars with M+ > 1.5 M©. 
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TABLE 1 

Properties of M Dwarfs with Giant Planets 



Gliese 


Hipparcos 


Spectral 


Stellar 


[Fe/H] 


M P sin i (M Jup ) 


Scmimajor 


Reference 


Number 


Number 


Type 


Mass (M ) 






Axis (AU) 




876 


113020 


M4 


0.32 


+0.37 


0.6189, 1.9275 a 


0.207, 0.130 


Rivera et al. (2005) 


849 


109388 


M3 


0.45 


+0.58 


0.83 b 


2.35 


Butler et al. (2006) 


317 




Ml 


0.24 




1.17 


0.95 


Johnson et al. (2007) 


832 


106440 


M2 


0.45 


-0.12 


0.64 


3.4 


Bailey et aT*T2009) 


179 


22627 


M4 


0.36 


+0.30 


0.9 


2.42 


Howard et al. (2009) 




79431 


M3.5 


0.50 


+0.4 


1.1 


0.34 


Apps et al. (2010, submitted) 


649 


83043 


Ml. 5 


0.54 


+0.1 


0.328 


1.135 


This work 



The G1876 planetary system contains two resonant Jovian planets and an inner "super-Earth" with Mp sini 

4.7 m s — 1 



b The orbit solution for G1849 includes a linear velocity trend dv/dt = 
planet. 

c The orbit solution for G1317 includes a linear velocity trend dv/dt 
planet. 



5.9 Af ffl . 



yr 1 , which may correspond to a second 
yr — , which may correspond to a second 



TABLE 2 
Stellar Properties and Orbital 
Solution for Gl649 



Parameter 


Value 


V 


9.70 ± 0.04 


K 


5.62 ± 0.02 


V -K 


4.08 ± 0.05 


B-V 


1.52 ± 0.04 


M v 


9.63 ± 0.05 


M K 


5.55 ± 0.02 


d (pc) 


10.34 ± 0.15 


M* (M s ) 


0.54 ± 0.05 


Teff (K) 


3700 ± 60 


[Fe/H] 


+0.08 ± 0.06 


P (days) 


598.3 ± 4.2 


P (years) 


1.638 ± 0.011 


K (ms^ 1 ) 


12.4 ± 1.1 


e 


0.30 ± 0.08 


T P (Julian Date-2400000) 


12876 ± 22 


u) (degrees) 


352 ± 15 


M P sin i (M Jup ) 


0.328 ± 0.032 


a (AU) 


1.135 ± 0.035 


N obs 


44 


rms (ms — 1 ) 


4.2 




1.17 




1.37 
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TABLE 3 
Radial Velocities for Gl 649 



JD RV Uncertainty 

-2440000 (m s" 1 ) (m s" 1 ) 



11409 


.824 


-1.54 


3. 


.35 


11705 


.913 


25.38 


3 


.80 


12004 


.048 


-8.26 


3 


73 


12007 


.010 


-10.05 


4. 


01 


12008 


.001 


-11.88 


3 


.92 


12009 


.070 


-2.72 


3 


65 


12097 


.968 


-8.16 


3 


.60 


12098 


.912 


1.04 


3 


.83 


12099 


.844 


4.90 


3 


63 


12100 


.895 


4.57 


3 


.39 


12127 


.911 


0.00 


3 


.72 


12161 


.808 


5.87 


3 


.89 


12189 


.748 


2.16 


4 


36 


12390 


.092 


13.12 


3 


,74 


12486 


.759 


-1.86 


3 


46 


12514 


.756 


-5.61 


3 


.56 


12535 


.722 


-6.90 


3 


.57 


12574 


.690 


2.29 


4. 


02 


12712 


.157 


-3.63 


3 


.71 


12777 


.039 


7.53 


3 


.93 


12804 


.920 


13.40 


3 


.75 


12832 


.961 


13.60 


3 


91 


13179 


.964 


-9.51 


3 


.82 


13195 


.812 


0.17 


3 


38 


13430 


.124 


13.32 


1 


20 


13547 


.904 


12.84 





.99 


13842 


.077 


-3.95 





97 


13934 


.875 


9.44 





95 


14640 


.044 


14.79 


1 


16 


14671 


.928 


16.81 





.72 


14673 


.897 


25.43 





.58 


14779 


.698 


8.75 


1 


,14 


14964 


.087 


4.24 


0. 


.96 


14985 


.917 


-3.91 





96 


15014 


.851 


-2.75 





.98 


15015 


.879 


-7.25 





.97 


15019 


.036 


-2.05 


1 


51 


15041 


.950 


-9.03 


1 


.05 


15042 


.786 


-6.49 





.97 


15043 


.900 


-9.16 


1 


00 


15048 


.857 


-3.67 


1 


.28 


15075 


.744 


4.40 





.91 


15111 


.736 


-2.10 


1 


.07 



